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We study the Reynolds number scaling of a gain-based, low-rank approximation to 
turbulent channel flows, determined by the resolvent formulation of McKeon & Sharma] 



(2010), in order to obtain a description of the streamwise turbulence intensity from 
direct consideration of the Navier-Stokes equations. Under this formulation, the veloc- 
ity field is decomposed into propagating waves (with single streamwise and spanwise 
wavelengths and wave speed) whose wall-normal shapes are determined from the prin- 
cipal singular function of the corresponding resolvent operator. We establish that the 
resolvent formulation admits three classes of wave parameters that induce universal be- 
havior with Reynolds number on the low-rank model, and which are consistent with 
scalings proposed throughout the wall turbulence literature. For the rank-1 model sub- 
ject to broadband forcing, the integrated streamwise energy density takes a universal form 
which is consistent with the dominant near- wall turbulent motions. When the shape of 
the forcing is optimized to enforce matching with results from direct numerical simula- 
tions at low turbulent Reynolds numbers, further similarity appears. Representation of 
these weight functions using similarity laws enables prediction of the Reynolds number 
and wall-normal variations of the streamwise energy intensity at high Reynolds numbers 
(Re T ~ 10 3 — 10 10 ). Results from this low rank model of the Navier-Stokes equations 
compare favorably with experimental results in the literature. 



1. Introduction 

Understanding the behavior of wall-bounded turbulent flows at high Reynolds numbers 
has tremendous technological implications, for example, in air and water transportation. 
This problem has received significant attention over the last two decades especially in 
the light of full-field flow information revealed by direct numerical simulations (DNS) 
at relatively small Reynolds numbers and high- Reynolds number experiments. Notwith- 
standing the recent developments, the highest Reynolds numbers that are considered in 
DNS are an order of magnitude smaller than experiments, which are in turn conducted 
at Reynolds numbers that are typically two orders of magnitude smaller than most ap- 
plications. This creates a critical demand for model-based approaches that describe and 
predict the behavior of turbulent flows at technologically relevant Reynolds numbers. 
Wall turbulence has been the topic of several revie ws; see, for example, [ Robinson 



(1991); Adrian (2007) for structure of coherent motions, Gad-El-Hak & Bandyopadhyay 
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(1994) for turbulence statistics and scaling issues, |Panton| (|2QQ1 ) for self-sustaining tur- 



bulence mechanisms, and |Klewicki| ( |2010| ); |Marusic et a/.| ( |2010c 



Smits et oZj ( [2011] ) for 



the latest findings and main challenges in examining high-Reynolds number wall turbu- 
lence. In the present study, special attention is paid to scaling and universality of the 
turbulent energy spectra at high Reynolds numbers. We also note that the energy spectra 
exhibit clear signatures of coherent turbulent motions such as the near- wall streaks, the 
large-scale motions (LSMs), and the very large-scale motions (VLSMs). 

1.1. Overview of dominant coherent motions 

In the interests of giving a brief overview of the energetically dominant coherent motions 
in wall turbulence, we will review three classes of structure. The near- wall system of quasi- 
streamwise streaks and counter-rotating vortices with streamwise length and spanwise 
spacing of approximately 1000 and 100 inner (viscous) units, centered at approximately 
15 inner units above the wall, has been well-studied. These ubiquitous features of wall 



turbulence are responsible for large production of turbulent kinetic energy ( [Kline et al. 
1967||Smith fc Metzlerp983| >. 



Another commonly-observed feature of turbulent flows is the hairpin vortex. In low 
Reynolds number flows, at least, packets of hairpin vortices have been observed to extend 
from the wall to the edge of the boundary layer and constitute LSMs (He ad fc Bandy- 1 
opadhyay 198l| [Adrian et aT]|2000 Adrian||2007 ), with streamwise extent approximately 
2-3 outer units (channel half-height, pipe radius, or the boundary layer thickness). 

VLSMs have been observed to reside in the logarithmic region of the turbulent mean 
velocity, with lengths of approximately 10-15 outer units in boundary layers and up to 
30 outer units in channels and pipes (see, for example, Kim fc Adrianp~999 Balakumar 



fc Adrian|[2007| |Monty et aZ.||2QQ7|). The emergence of VLSMs was originally attributed 



to alignment of LSMs ( Kim fc Adrian|[l999| ). However, |Smits et al\ ( |2011[ ) concluded 
that this is unlikely since the detached LSMs are located at a farther distance from the 
wall than the VLSMs and the attached LSMs have much smaller width than VLSMs 
and are convected at different speeds. Recently, the correlation between the envelope of 
small scale activity and the large scale velocity signal (identified via filtering in spectral 
space), which has been interpreted as an amplitude modulation of the small scales, has 
been investigated in detail, e.g. Hut chins fc Marusic] ( |20076| ); |Mathis et al. ( 2009a|6 ); 



Chung fc McKeon (2010); Hutchins et al (2011). 



1.2. Overview of scaling issues 

In spite of recent advances in understanding the structure of wall turbulence, the Reynolds 
number scaling of the turbulent energy spectra and the energy intensities remains an 
open area of research. The main experimental obstacle is maintaining the necessary spa- 
tial resolution for measurement accuracy while achieving the high Reynolds numbers 
required for large separation between the small and large turbulent scales. For example, 
the available experiments are performed at relatively small friction Reynolds numbers, 
Re T « ( 9(10 4 ), with a notable exce ption of the atmospheric surface layer measurements 



of, e.g., Metzger & Klewicki (2001) (Re T « O(10 6 )) that are in turn generally contami- 



nated by surface roughness effects. Most high-Reynolds number experiments suffer from 
spatial resolution issues in the inner region (see, for example, Hutchins et a/.|[2009[ ). 

Significant experimental effort has been devoted to determining the behavior of the 
streamwise energy intensity at high Reynolds numbers since it dominates the turbulent 
kinetic energy and is easier to measure relative to the wall-normal and spanwise velocities. 
It is understood that both small and large scales contribute to the streamwise energy 
intensity ( |Metzger fc Klewicki|200T| |Marusic fc Kunkel|2003[ |Hutchins fc Marusic|2007a 
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Marusic et al. 2010a). It is well-known that a region of the streamwise wavenumber 
spectrum scales with inner units; |Marusic et al.\ ( |2010a| ) showed by filtering that the 
contribution of such scales to the streamwise energy intensity, and therefore by extension 
also the streamwise spectrum, is universal, i.e. independent of Reynolds number. On the 



other hand, the large motions have been shown to scale in outer units (Kim & Adrian 



1999 ); Mathis et al. (2009 a ) proposed that the corresponding peak in streamwise intensity 



occurs close to the geometric mean of the limits of the logarithmic region in the turbulent 
mean velocity. The amplitude of this energetic peak increases with Reynolds number 
and has a footprint down to the wall (Hut chins & Marusic 20076). Note however, that 



available data are not sufficiently well-resolved to determine unequivocally the Reynolds 
number scaling of either the inner or outer peaks of the streamwise energy intensity (see, 
for example, Marusic et al. [20 10a ) 



Theoretical approaches also offer insight into the scaling of the spectrum with in- 
creasing Reynolds numbers, originating with the attached eddy concepts described by 



Townsend (1976). Perry & Chong (1982) developed these ideas to include hierarchies of 



geometrically similar, attached eddies in the logarithmic region, and formulated a model 
for the distribution of eddies that was required to capture characteristics including the 
mean velocity and streamwise turbulence intensity profiles. Subsequent works by Perry 
and co-authors extended this formulation beyond the wall region; most recently Ma rusic| 
|fc KunkeT] ([2003 ) used empirical scaling arguments concerning the effective forcing of the 
outer turbulence on the viscous region to propose a similarity expression for the stream- 
wise energy intensity that is valid throughout the zero pressure boundary layer. More 



recently, Marusic et al. (|20106| outlined an observationally-based, predictive formulation 
for the variation of the streamwise turbulent intensity up to the geometric mean of the 
logarithmic region based on consideration of the correlation between large and small 
scales. 

1.3. Review of previous model-based approaches 

We seek in this work a description of the streamwise turbulence intensity for all wall- 
normal locations arising from direct consideration, and modeling, of the Navier-Stokes 
equations. There has been much work in this vein, highlighting several important features 
of the Navier-Stokes equations. We provide a brief review of the most relevant literature 
here. 

It is understood that high sensitivity of the Navier-Stokes equations (NSE) linearized 
around the laminar velocity provides alternative paths to transition that bypass linear 
instability; see, for example, Schmid & Henningson (2001 ). Trefethen et al. ( |1993| ) showed 
that the high flow sensitivity is related to non-normality of the coupled Orr-Sommerfeld 
and Squire operators; see also Schmid ( 2007| ). These operators are coupled in the pres- 
ence of mean shear and spanwise- varying fluctuations. Physically, as originally explained 
by Landahl (19751), a large streamwise disturbance is induced on the flow in response to 



lift-up of a fluid particle by the wall-normal velocity such that its wall-parallel momentum 
is conserved. 

Even in linearly stable flows, the high sensitivity can result in large transient responses, 
meaning that the energy of certain initial perturbations significantly grows before even- 



tual decay to zero (jGustavsson 1991 Klingmann 1992 Butler & Farrell 1992 Schmid 



fc Henningson|[l992| Reddy fc Henningson||1993 ). In addition, the high sensitivity is re- 



sponsible for high energy amplification, meaning that the velocity fluctuations achieve 
a large variance at the steady state for the flow subject to zero-mean stochastic distur- 



bances (Farrell & Ioannou 19936 Bamieh & Dahleh 2001 Jovanovic & Bamieh 2005). 



The dominant structures that emerge from the above transient growth and energy am- 
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plification analyses are reminiscent of the streamwise streaks observed at the early stages 
of transition to turbulence ( jMatsubara fc Alfredsson]|2QQl| ). They are characterized by 
infinitely long spanwise-periodic regions of high and low streamwise velocity surrounded 
by pairs of counter-rotating streamwise vortices that are separated by approximately 3.5 
outer units. 

The critical role of linear amplification mechanism in promoting and maintaining tur- 
bulent flows was highlighted in direct numerical simulations of Kim & Lim (2000). It 
is known that the NSE linearized around the turbulent mean velocity are stable for all 
Reynolds numbers (Malkus 1956). Early model-based approaches extended the afore- 
mentioned sensitivity analyses of the laminar flow to the turbulent channel flow and 
found dominance of streamwise streaks that are spaced by 3 outer units, which is ap- 
proximately the same as in the laminar flow. In addition to the outer-scaled dominant 
structures, |Butler fc Farrell (1993); Farrell & Ioannou (1993a]) showed that the largest 
transient response over an eddy turnover time of 80 inner units, associated with the 
near-wall cycle, is obtained for initial perturbations that are infinitely long and have the 
same spanwise spacing as the near-wall streaks, i.e. 100 inner units. The same stream- 
wise and spanwise lengths were obtained in flows subject to stochastic disturbances over 
a coherence time of 90 inner units ( Farrell fc Ioanno"u| [l998). 

Reynolds fc Hussain| ( |1972[ ) put forward a modified linear model to account for the 



effect of background Reynolds stresses on the velocity fluctuations. They proposed to aug- 
ment the molecular viscosity by the turbulent eddy viscosity that is required to maintain 
the mean velocit y. This model yields two local optima for the structures with largest 



transient growth ( del Alamo fc Jimenez||2006 Pujals et a/.| |2009) and energy amplifica- 
tion ( Hwang fc Cossu|2010 ) without the need for confining the optimization time. These 
peaks correspond to streamwise-elongated structures with a spacing of 80 inner units and 
3-4 outer units and are in fair agreement with the spacing of near-wall streaks and the 
very large-scale motions in real turbulent flows. 

An exact representation of the NSE was introduced by McKeon fc Sharma] (2010) in 
which (i) a set of linear sub-systems describe extraction of energy from the mean velocity 
at individual wavenumbers /frequencies; and (ii) the only source of coupling between these 
sub-systems is the conservative nonlinear interaction of their outputs, that determines 
both the input to the sub-systems and the turbulent mean velocity. At its heart is the 
ability to analyze the flow of energy from the mean velocity to all the velocity scales and 
identify the essential linear amplification and nonlinear redistribution mechanisms that 
drive the turbulent flow. The input-output relationship of the linear sub-systems can 
be described by transfer functions whose low-rank nature in the wall-normal direction 
enables significant simplification of their analysis. 



One of the main differences between the formulation of McKeon & Sharma (2010) and 



other input-output analyses of laminar and turbulent flows (see, for example, 



Jovanovic 



fc Bamieh|2005 Hwang fc Cossu|2010 ) is parameterization of the waves with wave speed 



rather than temporal frequency. While the latter approaches conclude that the glob- 
ally optimal transient growth and energy amplification takes place for zero streamwise 
wavenumber and temporal frequency, they do not provide guidelines for finding locally 
(in wall-normal direction) optimal wave parameters. On the other hand, selecting the 
wave speed (i) enables systematical search for both locally and globally optimal wave 
shapes and parameters; (ii) removes the ambiguity about the wave speed corresponding 
to the globally optimal waves by determining the limit of the ratio between zero stream- 
wise wavenumber and temporal frequency; and (iii) distinguishes between non-normality 
and criticality as the main linear amplification mechanisms. 



McKeon & Sharma ( 2010[ ) showed that the principal forcing and response directions 
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associated with the linear sub-systems are consistent with the dominant response shapes 
in real turbulent pipe flows. In addition, the low-dimensional and sparse feature of the 
resulting model enables development and utilization of compressive sampling techniques 



for analyzing the turbulent flow dynamics ( |Bourguignon et al. 2013). This formulation 
has also proven useful for pre- and post-diction of experimental observations in turbulent 
pipe flow QSharma fc McKeon|[2Q12l |McKeon et a/.||2Q13 ). 



1.4. Paper outline 

In this paper, we identify the Reynolds number scaling of a low-rank approximation to 
turbulent channel flow and utilize it for predicting the streamwise energy intensity at 
high Reynolds numbers. Our development is outlined as follows: In §[2j we briefly review 
the resolvent formulation, highlight its low-rank nature, and show that a rank-1 approx- 
imation captures the most energetic modes of real turbulent channels. The stage is set 
for studying the energy density of fluctuations using a minimum number of assumptions 
by considering a rank-1 model in the wall- normal direction subject to broadband forcing 
in the wall-parallel directions and time. Furthermore, a summary of the computational 
approach for determining the rank-1 model is provided. 

Three classes of wave parameters for which the low-rank approximation of the resolvent 
exhibits universal behavior (independence) with Reynolds number are identified in § [3] 
The requirement for universality highlights the role of wave speed in distinguishing these 
classes. Each class of waves is characterized by a unique range of wave speeds and a 
unique spatial scaling that emerge from the resolvent. For the rank-1 model subject to 
broadband forcing, we reveal the universal streamwise energy densities, and show that the 
peaks of these energy densities roughly agree with the most energetic turbulent motions, 
i.e. the near-wall streaks, the VLSMs, and the LSMs. 

In §|4j we show that the streamwise energy density of the rank-1 model with broadband 
forcing can be optimally weighted as a function of wave speed to match the intensity of 
simulations at low turbulent Reynolds numbers. The weight functions are then formulated 
using similarity laws which, in conjunction with the universal energy densities, enable 
prediction of the streamwise energy intensity at high Reynolds numbers. The paper is 
concluded in § [5] and several future directions are discussed. 



2. Low-rank approximation to channel flow 

An overview of the rationale for considering a low rank approximation to turbulent 



channel flow is presented in this section. We follow the development of|McKeon fc Sharma 



(2010) for turbulent pipe flow, showing that equivalent results are obtained for channels 
and highlighting the new observations. The pressure-driven flow of an incompressible 
Newtonian fluid is governed by the nondimensional NSE and the continuity constraint 

u t + (u • V)u + VP = (l/#e T )Au, 

Vu = 0, (2 - 1} 

where u(x, y, z, t) is the velocity vector, P(x, y, z, t) is the pressure, V is the gradient 
operator, and A = V • V is the Laplacian. The streamwise and spanwise directions, x 
and 2, are infinitely long, the wall-normal direction is finite, < y < 2, and t denotes 
time; see figure [T(aJ| for the geometry. The subscript t represents temporal derivatives, 
e.g. ii£ = du/dt. The Reynolds number Re T = u T hjv is defined based on the channel 
half-height h, kinematic viscosity z/, and friction velocity u T = \/r w /p, where r w is the 
shear stress at the wall, and p is the density. Velocity is normalized by u T , spatial variables 
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(b) 

Figure 1. (a) Pressure driven channel flow; and (b) Schematic of a two-dimensional propagating 
wave with streamwise and spanwise wavelengths X x and X z and streamwise propagating speed 
c. 

by h, time by h/u T , and pressure by pv%.. The spatial variables are denoted by + when 
normalized by the viscous length scale v/u T , e.g. y + = Re T y. 

2.1. Decomposition in homogenous directions 

The velocity is decomposed using the Fourier transform in the homogenous directions 
and time 



u(x,y,z,t) = JJJ u(y]K x ,K z ,uj)e [ ^ xX + KzZ ^d^d^dcj, 



(2.2) 



where A denotes a variable in the transformed domain, and the triplet (k x ,k, Z: uj) is the 
streamwise and spanwise wavenumbers and the temporal (angular) frequency. For any 
(k x , n z1 uo) 7^ 0, u(y; k x , k z ,uj) represents a propagating wave with streamwise and span- 
wise wavelengths \ x — 2tt / ' n x and \ z = 2tt / ' k z and speed c = uo/n x in the streamwise 
direction; see figure l(b)| for an illustration. Some special cases include standing waves 
(c = 0), infinitely long waves (k x = 0), and infinitely wide waves (k z =0). In this study, 
we emphasize the eminent role of wave speed, a factor that was highlighted by |McK-| 
|eon fc Sharmaj ( |2Q1Q| ) while being predominantly neglected in the previous studies, in 
determining the classes of propagating waves that are universal with Reynolds number. 

The turbulent mean velocity U(y) = [U(y) 0] T = u(y; 0,0,0) corresponds to 
(k, x ,k z ,uj) = and is assumed to be known. |McKeon fc Sharma ( 2010| ) avoided the 



closure problem for the mean velocity by using U (y) obtained in pipe flow experiments, 
but note that the resolvent formulation could be used to determine the mean velocity 
profile, a topic of ongoing work (see |McKeon et al. 2013). Here, we use a semi-empirical 
turbulent viscosity model, originally proposed for pipe flow (Malkus 1956| |Cess |1958 ) 
and extended to channel flow (Rey nolds fc Tiederman|[l967 ), to determine U(y) 



U(y) 



vr{y) 




4y + 2y 2 ) { 1 



(|y-i| 




1 

" 2' 
(2.3) 

where vt is normalized by ^, and the parameters a and n appear in the van Driest's 
wall law and the von Karman log law. These parameters are obtained by minimizing 



the deviation between U (y) in (|2.3|) and the DNS-based turbulent mean velocity profile. 
The a and k obtained by Moarref & Jovanovic ( 2012[ ) suggest that both of these values 
converge for large Re T . We take a = 25.4 and k = 0.426 for all Reynolds numbers and 
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note that these values are optimized for Re T = 2003 jdel Alamo fc Jimenez||2006 Pujals 
et fl/.|2009[ ). 



Following McKeon & Sharma (2010), the convective nonlinearity in (2.1) is considered 



as a forcing term f = — (u- V)u that drives the velocity fluctuations. For any (k x ,k z ,w) ^ 
0, an equation for velocity fluctuations u(y; k x , k z ^uj) = [u v w] T around the turbulent 
mean velocity is obtained by substituting (2.2) in ( |2.1| ), and using the orthonormality of 
the complex exponential functions 



-iuju + (U • V)u + (u • V)U + Vp - (l/i?e r )Au 
V-u = 0. 



(2.4) 



Here, V = [w x d y m z ] T , A = d yy — n 2 with n 2 = k x + k 2 : and 

f(y;K x ,K z ,uj) = [fi f 2 h] T = 

(u(y; k' x , k' z ,w') • V) u(y; k x — n f x , k z — n f z , uj — u/) &k! x dnc f z dou' '. 



(4,4,^)^(0,0,0) 

(4,4,u/) + (Kx,Kz,u) 

(2.5) 

McKeon & Sharma (2010|) implicitly accounted for the continuity constraint by pro- 



jecting the velocity field onto the divergence-free basis of |Meseguer fc Trefethen (2003). 



Here, we use a standard choice of wall-normal velocity v and wall-normal vorticity 
fj = m z u — 'm x w as the state variables, £(?/; k x , k Z: lj) = [v t)] t , to eliminate the pressure 
term and the continuity constraint from (2.4) and obtain 



-(lUjI + A(K, x ,K>z))C(y]Kx,K>z,v) = Ci(K, x ,K, z )f(y]K, x ,K, z ,U)), 
il(y;K x ,K z ,Lj) = C(K, x ,K, z )C(y]Kx,K z ,Cj). 



(2.6) 



Here, A is the state operator, C maps the state vector to the velocity vector, and the 
adjoint of C (denoted by C^) maps the forcing vector to the state vector. A, £?, and C 



are operators in y and parameterized by k x and k z \ see (3.1), to be discussed later. The 
input-output relationship between f and u is obtained upon elimination of £ from (2.6) 



u(y;n x ,K z ,uj) = H(k x , k z , uj) f (y; k x , k z , u), 



(2.7) 



where Ra(k x , k z ,u) = —(\ujI + A(k x , k z ))~ 1 is the resolvent of A. As illustrated in fig- 
ure^ the only source of coupling between propagating waves with different wavenumbers 
is the quadratic dependence of {(x,y,z,i) on u(x,y,z,i). For any wavenumber triplet, 
the input-output map from f to u (shown by the dashed rectangle) represents a sub-unit 
of the full NSE. 



2.2. Decomposition in the wall-normal direction 

The transfer function H{n xi n z1 uo) provides a large amount of informatio n about the 
input-output relationship between f and u. Following the gain analysis of 



McKeon & 



Sharma (2010), we use the Schmidt (singular value) decomposition to provide a wall- 
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FT 



f = _( u . V)u 
nonlinearity => coupling 



IFT 



Figure 2. For any triplet (k x , k z , cj), the operator H(k x , k z , uj) maps the forcing f to the 
response u. The different wavenumbers are coupled via the quadratic relationship between 
f (x, y, z, t) and u(x, y, z, t). FT and IFT stand for Fourier transform and inverse Fourier trans- 
form, respectively. The input-output map (shown with the dashed rectangle) is the main focus 
of the present study 



normal basis based on the most highly amplified forcing and response directions: 

CO 



(2.8) 



where ai > <r 2 > • • • > denote the singular values of H, and the singular functions 
4>j — [fij hj hj] T an d ^Pj — [% % ^j] T are respectively the forcing and response 
directions corresponding to aj. 

While the singular values of H are unique, additional treatment is necessary to obtain 
unique singular functions. Unlike in a pipe, the singular values come in pairs due to the 
wall-normal symmetry in the channel; see, for example, figure [4(a) The paired singular 



values are equal unless both the streamwise and spanwise wavelengths are large relative 
to the channel height (0(X X ) > 1 and 0(X Z ) > 1). When the paired singular values 
are different, unique singular functions can be obtained modulo a complex multiplicative 
constant of unit magnitude if an orthonormality constraint is imposed on them 



J ^]{y]n x ,K z ^)^ k (y\n x ,n z ,uj)dy = J 



^ j (y;K JXl K z ,uj)tl) k (y' 1 K JXl K JZ ,uj)dy = S jk , 



(2.9) 

where S denotes the Kronecker delta. In the case where the paired singular values are 
equal, we impose a symmetry /anti-symmetry constraint on the singular functions in ad- 
dition to the above orthonormality constraint. In other words, the corresponding singular 
functions assume the same magnitude throughout the channel while being in phase in 
one half of the channel and out of phase in the other half. 

In this study, we select the unknown multiplicative constant (after orthonormalization) 
such that Uj(y m &x', k x , k z ,ou) is a real number at the wall-normal location y max where the 
absolute value of Uj is the largest. This choice places the maximum of Uj(x : y, z,t; k x , k Z: oj) 
at the origin x = z = t = 0. The channel symmetries in the streamwise and spanwise 
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Figure 3. The principal velocity response i/) 1 (x,y, z,t; k x , k z ,oo) — [ui v\ w\] T for A+ = 700, 
Aj = 100, c = 10, and Re T — 10000 at t = 0. (Left) The isosurfaces of the streamwise velocity 
at 60% of its maximum; (Right) The streamwise velocity (contours) and the spanwise and 
wall- normal velocity (arrows) at x + — Aj/2. The contours in (b) represent positive (thick solid) 
and negative (thin dashed) values from 3 to 15 with increments of 3. 



directions can be used to obtain Uj, Vj, and Wj in the physical domain 

Uj(x,y, z,t; k x , k Zj uj) = 4 cos(k z z) Re (uj(y; k x -> k>z, w) e 1 ^ 3 ^ - ^) 

Vj(x,y,Z,t]K, x ,K z ,Uj) 
Wj(x,y,Z,t]K x ,K, z ,U)) 

where Re and Im denote the real and imaginary parts of a complex number. The rep- 
resentation of the forcing directions in the physical domain is obtained using similar 
expressions. 



4cos(^)Re (v j (y]K x ,K, z ,u)e i ( KaiX , 
— 4sin(K z z) Im (wj(y; k xj k z ,uj) Q^xx-ut)^ 



From the singular value decomposition (2.8) and the orthonormality constraints (2.9) 



it follows that if the forcing is aligned in the ■ direction with unit energy, the response 
is aligned in the ipj direction with energy <j|. Consequently, the forcing and response 
directions with the largest gain correspond to the principal singular functions <p 1 and 
jp 1 . For any (k, x , k z ,lj), the singular functions of H should be thought of as propagating 
waves in the physical domain. 

Equivalent near- wall structures to those reported for pipe flows by |McKeon fc Sharma 
( 2010t ; |McKeon et al.| ( |20l3| are obtained for channel flows. For example, the principal 



singular function ipi(x, y, z, t; k x ^ k z , u) = [u\ v\ w\] T for the propagating wave corre- 
sponding to the energetic near-wall cycle (A+ = 700, A+ = 100, c = U(y + = 15) = 10) 
for Re r = 10000 is shown in figure [3] The streamwise component of these structures 
contains regions of fast- and slow- moving fluids that are aligned in the streamwise direc- 
tion, slightly inclined to the wall, and are sandwiched between counter-rotating vortical 
motions in the cross-stream plane. 

2.3. Low-rank nature of H 

The operator H, acting on functions of y, can be described as low-rank if a significant 
portion of its response to a broadband forcing in y is captured by projection on the first 



few response directions. |McKeon fc Sharma (2010) highlighted the low-rank nature of 
H for turbulent pipe flow. Figure 4(a) shows the first twenty singular values of H for 
A+ = 700, A+ = 100, and c = 10 in turbulent channel flow with Re T = 2003. We see that 
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Figure 4. (a) The twenty largest singular values of H for A+ = 700, A+ = 100, c = 10, and 
Re T — 2003. (b)-(d): The color plots show the energy that is contained in the largest two 
response modes relative to the total response for different streamwise and spanwise wavelengths 
and (b) c = U(y + = 15); (c) c = U(y + = 100); and (d) c = U(y 0.2). The black c ontours 
show the turbulent kinetic energy spectrum from the DNS of Hoy as &; Jimenez] ( |2006| at the 
corresponding critical wall- normal locations (b) y + — 15; (c) y + — 100; and (d) y — 0.2. The 
contours represent 10% to 90% of the maximum energy spectrum at each wall-normal location 
with increments of 20%. 



the largest pair of singular values is approximately one order of magnitude larger than 
the other singular values. 

The energetic contribution of the fc-th direction ip k to the total response in the model 
subject to broadband forcing in y with fixed \ x , \ Z1 and c is quantified by a\j (X^j=i a j)- 
Figures 4(b)|4(d)| highlight the low-rank nature of H by showing that the first two princi- 
pal response directions ip 1 and ip 2 contribute to more than 80% of the total response over 
a large range of wall-parallel wavelengths (red region) for wave speeds c = U(y + = 15), 
U(y + = 100), U(y = 0.2), and Re T = 2003. The relevance of studying the low-rank 
approximation of H is further emphasized by noting that the wavenumbers for which 
H is low-rank coincide with the most energetic wavenumbers from the DNS of Hoy as fc 



Jimenez (2006) (contours) at the corresponding critical wall-normal locations. 
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2.4. Rank-1 model subject to broadband forcing 

In the present study, we consider a rank-1 model by only keeping the most energetic 
forcing and response directions corresponding to g\ and show that significant under- 
standing of the scaling of wall turbulence can be obtained using this simple model. This 



is motivated by the observation in § 2.3 that the operator H is essentially a directional 
amplifier. In other words, we expect to see the principal singular response of H in real 
turbulent flows provided that the principal forcing direction is present in the nonlinear 
forcing term; i.e. the gradient of the Reynolds-stresses. In order to use the least number 
of assumptions, we consider the case where the forcing f equals the principal forcing 
direction 1 . Consequently, the forcing has unit energy for all wave parameters, meaning 
that it is broadband m tz x i i 

and c. 

It is well-known that the streamwise energy spectrum can be divided into regions that 
scale in inner and outer variables (see, for example, Morrison et a/.|2004 ). Our objective is 



to explore the Reynolds number scaling of the streamwise energy density and predict the 
behavior of the streamwise turbulence intensity at high Re T . We focus on the streamwise 
velocity because it dominates the kinetic energy density in turbulent flows. Similarly, 
the principal singular responses of H that result in the largest energy amplification are 
dominated by their streamwise component. For the rank-1 model with broadband forcing, 
we define the premultiplied streamwise energy density of the principal response of H by 

Euu{W,Kx,K>z,c) = K 2 X K Z ((Ti(K x ,K z ,c)\ui\(y, K x ,K z ,C)f , (2.10) 

such that the premultiplied one-dimensional energy densities and the energy intensity 
are obtained by integrating E uu (y; k x , k z ,c) over the set of all wave parameters <S, e.g. 



E uu (y,c) = / / E uu (y; c)dlog(>c x )dlog(K 2 ), 

(2.11) 



s 

E uu (y) = J J J E uu (y-K x , K z ,c)dlog(K x )dlog(K z )dc, 
s 

and E uu (y : K x ) and E uu (y,K z ) are determined similarly. 

The above formulation of the energy density is used in § [3] to identify the contribution 
of confined subsets of wave parameters to the energy density. We establish that the 
energy density exhibits universal behavior with Re T for properly selected subsets of 
wave parameters. It is further shown that the emerging scales are consistent with those 
observed in experiments. In addition, the scales of energetically dominant waves roughly 
agree with the scales of dominant near-wall motions in real turbulent flows. 

2.5. Computational approach 

A pseudo-spectral method is used to discretize the differential operators in the wall- 
normal direction on a set of Chebyshev collocation points. This is implemented using the 
Matlab Differentiation Matrix Suite developed by Weideman & Reddy (2000). Table [I] 



summarizes the selected range of wave parameters and their respective resolution in 
numerical computations. It has been verified that the excluded wave parameters are not 
energetically important and therefore do not change the results of the present study. 



An efficient randomized scheme developed by Halko et al. (2011) is utilized to compute 
the principal singular directions of H for different Reynolds numbers and wave parame- 
ters. The accuracy and computation time depend on the decay of the singular values; a 
faster decay results in high accuracy or equivalently less computation time to reach the 
same accuracy. In addition, if the singular values are not well separated, the problem 
of computing the associated singular functions is badly conditioned, meaning that it is 
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Re T 


N x 


Ny 


N z 


N c 


\ + 

cc,min 


Acc,max 


y + . 


\ + 

\z,min 


A^max 


U c i 


934 


64 


251 


32 


100 


10 


10 6 


0.07 


10 


100 


22.39 


2003 


64 


251 


32 


100 


10 


5 x 10 5 


0.15 


10 


50 


24.02 


3333 


64 


301 


32 


100 


10 


3 x 10 5 


0.18 


10 


30 


25.22 


10000 


64 


401 


32 


100 


10 


10 5 


0.3 


10 


10 


27.81 


30000 


80 


601 


40 


100 


10 


3.3 x 10 6 


0.4 


10 


33 


30.39 



Table 1 . Summary of the selected parameters in numerical computations at different Reynolds 
numbers. In the wall-normal direction, N y Chebyshev collocation points is used with y^ in de- 
noting the closest point to the wall. In the streamwise and spanwise directions, N x and N z 
logarithmically spaced wavelengths are used between A+ in and A max . In addition, N c linearly 
spaced wave speeds are chosen between c m i n = 5 and c max = U c i> 



hard for any method to determine them very accurately. In this study, the above scheme 
approximately halves the total computation time relative to Matlab's svds algorithm. 
This becomes increasingly important considering the three-dimensional wave parameter 
space that we need to explore and the large size of the discretized resolvent operator 
(twice the number of collocation points in y) at high Reynolds numbers. In addition, 
the randomized nature of this scheme enables its parallel implementation which makes 
it especially suitable for large-scale computations. Even though we have not used this 
feature in the present study, it may find use in designing turbulent flow control strategies, 
e.g. by means of spatially or temporally periodic actuations. 



3. Reynolds number scaling 

The formulation of § [2] facilitates analysis of the contribution of different wave param- 
eters (k x ,k z ,c) to the streamwise energy density. For the rank-1 model with broadband 
forcing, the energy density of each wave is determined from the principal singular values 
and singular functions of the transfer function H; see ( 2.10| ). In this section, we identify 



three classes of wave parameters for which E uu (y; k x , k z ,c) exhibits universal behavior 
with Re T . Each class is characterized by a unique range of wave speeds and a unique 
scaling of the wall-parallel wavelengths. 



3.1. Requirement for universal behavior 

We start by showing that a requirement for universal behavior is the wall-normal locality 
of the singular functions of H. This is done by examining the underlying operators in H, 



cf. (2.7): 



A = 
C = 
Ra = 



1 



A" 1 ((l/Re T )A 2 + m x (U" - UA)) 

—m z U' (l/Re T )A - m x U 



K 

W Z dy 



ik x A 1 d y k 2 A 1 -i^A x dy 



\K Z 







A" 1 (ik x (([/ - c)A - U") - (l/i?e T )A 2 ) 

\k z U' \k x {U-c) - (l/i?e r )A 



(3.1) 
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a, 



yyyy 



yy 



and prime denotes differentiation in y, e.g. U'(y) 



where A 2 

dU/dy. We see that the difference between the turbulent mean velocity and the wave 
speed, U(y) — c, and its wall- normal derivatives, U'(y) and U"(y), appear as spatially- 
varying coefficients in H. Since the turbulent mean velocity scales differently with Re T 
in different wall-normal locations, only the singular functions that are sufficiently narrow 
in y have the potential to be universal. This is because such singular functions are purely 
affected by a certain part of the mean velocity that scales uniquely with Re T . 

We next show that the singular functions c orresponding to the energet ically signifi- 
cant modes are in fact localized. Recent DNS ( del Alamo fc Jimenez||2009 ) and experi- 
ments ( LeHew et a/.|201l| ) have shown that the energetic contribution of structures with 
convection velocities less than 10u T and larger than the centerline velocity U c \ = U(y = 1) 
is negligible in real turbulent flows. This motivates defining a conservative subset of <S, de- 
noted by <S e , that includes all wall-parallel wavenumbers and the energetically important 
wave speeds 

S e = {(k x ,k z ,c) I 5<c<U c i}. (3.2) 

Since a significant portion of the streamwise energy intensity is contained below y + ~ 15 
(critical speed c « 10), we include the waves with c < 10 in S e such that the principal 
singular functions peak for y + < 15. 

Figure [5] shows the one-dimensional energy density as a function of wave speed E uu (y, c) 
for Re T = 30000 and S = S e . As evident from figure [5(a) , the energy density for a fixed c 
is localized in a narrow wall-normal region; note that the colors are given in logarithmic 
scale. The localization is highlighted in figure [5(b)] where E uu (y,c) is normalized by its 
maximum value over y for fixed values of c. We see that the largest energy amplification 
takes place in the vicinity of the critical wall-normal location where the turbulent mean 



velocity (thick white curve) equals the wave speed. McKeon & Sharma (2010) argued 



that emergence of critical layers is one of the three means of maximizing the Hilbert- 
Schmidt norm of H (sum of squares of the singular values), i.e. by locally minimizing the 
term U(y) — c that appears in the resolvent operator Ra given in (3.1). 

According to Taylor's frozen turbulence hypothesis ( Taylor|[l938 ), the flow structures 
in boundary layers propagate downstream with a speed close to the local mean velocity. 
Consistent with this hypothesis, figure [5(b)] shows that among all the waves with arbitrary 
streamwise and spanwise wavelengths at the wall- normal location y, the ones with critical 
speed c = U(y) are the most highly amplified. This provides strong evidence for the 
importance of critical layers in amplification of flow disturbances. In addition, figure [5(b)] 
shows that the scatter in the energetic wave speeds increases as the peak of energy density 
approaches the wall. This agrees with the practical observation that Taylor's hypothesis 



yields inaccurate energy spectra close to the wal l : see, f o r example , Kim & Hussain ( 1993[ ); 
Monty fc Chong|(|2009|);|del Alamo fc Jimenezj (2009|); |LeHew et al\ $0u\. 



3.2. Summary of the universal behavior with Re T 

We start by reviewing the universal behavior of the turbulent mean velocity. This is a 
prerequisite to studying the universality of the principal propagating waves since the 
latter holds for critical modes only, as discussed in § |3.1| In the commonly accepted 



picture (Coles 1956), the mean velocity is divided into inner, logarithmic, and outer 



regions 



U = B(y+) + (1/k) log(y+) + (2U/ K )W(y), 



(3.3) 



where B is the inner-scaled wall function, II is the wake factor, W is the outer-scaled 
wake function, and the von Karman's constant n = 0.426 is taken to match the mean 
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Figure 5. (a) The one-dimensional energy density Re~ 2 E uu (y, c) for S = S e and Re T = 30000; 
and (b) the energy density normalized by its maximum value over all y for fixed values of c. The 
colors are in logarithmic scale. The turbulent mean velocity is shown by the thick white curve 
in(b). 



Region 
Inner 



Range of wave speeds 
c < 16 



Logarithmic \U m - c\ < log (-^ 6 ,mm 



50 



Outer 



2k "° v 50/0.2 
< Ud - c < 6 



Wall- normal scale 

y + < ioo 
y > o.i 



< 



Table 2. Summary of wave speeds and wall-normal intervals where U — c is universal with Re T . 
The mean velocity is denoted by U c i in the center plane and by U m in the geometric mean of 
the logarithmic region. The smallest Reynolds number that is considered is denoted by i?e r , m in- 
See also the green regions in figures 6(b)|6(d) 



Class 
Inner 


x-scale 


y-scale 

y + 


z-scale 
A+ 


t-scale 

- 1 C 


Si 


Outer-in-x 


Ax 






T 




Outer-in-2/ 


A x 
Re T 


y 


A, 


Re T T c i 


S oy 



Subset of wave parameters 
5 < c < 16 

T,min J 

< Ud - c < 6 

X x /X z > Vld(Re T /Re r,min ) 



Table 3. Summary of universal wave parameters with Re T . See also figures 6(e)|6(g) 



velocity given in (2.3) for high Re T . Consequently, U — c is universal with Re T for certain 
intervals of wave speed and wall-normal scales; see tableland figures 6(a)|6(d)| 

Figure 6(b) shows that U(y + ) — c is universal for y + < 100 and fixed c < 16 (inner 
region). As shown in figure 6(d) , the function U(y) — c is universal for y > 0.1 and fixed 
defect wave speeds relative to the centerline < U c \ — c < 6, with U c \ = U(y = 1) (outer 
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( outer-in-y 

\ x /Re T ,y,\ z 


S oy 








outer-in-x 


Sox 






^, \fyy+, vXaJ 








inner 


Si 


INNER 


LOG 


OUTER 





100 



0.1 




Figure 6. (a) Schematic of the the different regions of the mean velocity and the associated 
classes of Reynolds number scales that are induced on the propagating waves, (b) The turbulent 
mean velocity U(y + ); (c) the defect velocity U c i — U(y) defined relative to the centerline; and 



(d) the defect velocity U ( 



Um defined relative to the geometric mean of the logarithmic 
region y m y^ = 10, for Re T = 3333 (blue), Re T = 10000 (red), and Re T = 30000 (black). 
Notice that the shaded regions are invariant with Re T . (e)-(g): The classes of wave parameters 
(Ax, \ z , c) that induce universal behavior with Re T on the low-rank approximation of H. (e) Si 
induces the inner scale; (f) S oy induces the outer-in-y scale; and (g) S ox induces the outer- in-x 



scale. The classes S oy and S ox are obtained for Re T > Re T , 



3333. 



region). The logarithmic region, take n for y + > 50 and y < 0.2, bridges the gap between 
the inner and outer regions. Figure 6(c) illustrates that the function U {\fyy + ) — c is 
universal in the logarithmic region associated with the smallest Reynolds number Re TjYn i n 
that is considered, and for fixed defect wave speeds relative to the geometric mean of 



the logarithmic region \U n 
VTO). For example, Re TyYn 



■c\ < (1/2/c) logCRe^ 
= 3333 yields \U m — 



./(50/0.2)), with U m = U(, 



< 3 and 1 < 



10, i.e. y is 
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Class 


Subset 




(d/dy) 






Pi 


Ui 


E uu 


Inner 


Si 


Re T 


Re T 


Re T 


Re T 


Re' 1 


ReV 2 


Re 2 T 


Outer-in-x 


Sox 


1 


ReV 2 


ReV 2 


1 


ReV 2 


ReV 4 


Re 2 T 


Outer-in-y 


Soy 


Re' 1 


1 


1 


Re' 1 


Re 2 T 


1 


Re 2 T 



Table 4. Summary of the growth/decay rates with Re T of the wall-parallel wavenumbers, the 
wall- normal derivative, the convective frequency, the principal singular value and the principal 
streamwise singular function of if, and the premultiplied three-dimensional streamwise energy 
density for the three classes of universal waves with Re T . 



approximately within a factor of v 10 from the geometric mean of the logarithmic region. 
We note that the latter scale in the logarithmic region also emerges from the resolvent 
operator such that an outer scaling is imposed on the streamwise wavelength. 

Associated with each region of the mean velocity, there is a class of wave parameters 
for which the low-rank approximation of H exhibits universal behavior with Re r ; see 
tables [3] and [4] for a summary. As illustrated in figures [6(a)] and [6(e)|6 (g) | these classes are 
primarily distinguished by the wave speed. The principal singular functions inherit their 
wall- normal length scale from the turbulent mean velocity at the critical layer, and their 
streamwise and spanwise length scales from the resolvent. In addition, the magnitude 
of the singular values and singular functions scale uniquely with Re T in each class of 
wave parameters, which induces unique scales on the premultiplied streamwise energy 
density E uu (y; k X: k z , c). Next, we separately discuss each class and refer the reader to 
Appendix [A] for detailed derivation of the scales. 

3.2.1. The inner class Si 

For wave speeds in the inner region of the turbulent mean velocity, universality of H 
requires constant A+, ?/ + , A+, and c; cf. tableland figure 6(e) As a result, the time 
T c = A^/c over which the wave convects downstream for one wavelength relative to the 
wall reduces with Re T and the convective frequency uj c = 2tt/T c increases with Re T . In 
other words, a truly inner scale is induced on the length, height, width, and convective 
time of the waves that correspond to the principal singular functions. Therefore, the wall- 
normal support of the singular functions in outer units linearly decreases with Re T1 and 
the unit energy constraint on the singular functions requires that the magnitude of the 
singular functions increase with Re}/ 2 . The number of these waves per unit wall-parallel 
area and time increases with Re^ T as their length, width, and convective time decrease 
with Re T . Since the singular values of H linearly decrease with Re T1 the overall result is 
that E uu (y; k x , k, z ,c) increases with Re 2 .; cf. table [4] 

3.2.2. The outer-in-x class S ox 

For wave speeds in the logarithmic region of the turbulent mean velocity, universality 
of H requires constant A^, y + y, A+A 2 , and U m — c, such that the a spect ratio constraint 
^x/^z ^ \/^0Re T /Re TiYn i n is satisfied; cf. tableland figure 6(f) As a result, the time 
T m = ^x/\U m — c\ over which the wave convects away for one wavelength relative to 
an observer with speed U m remains independent of Re T and same does the convective 
frequency uj m = 2ir /T m . These waves have the same scales as the structures in the meso- 



layer; see, for example, |Long fc Che n (1981); Afzal ( |1984 ); Sreenivasan & Sahay (1997); 
Wei et al. ( 2QQ5| ). We see that the aspect ratio \ x /\ z of the universal waves increases 
as This explains why universality for this class holds for the waves with aspect 
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ratios larger than a threshold: As the aspect ratio of the singular modes increases with 
Reynolds number, the Laplacian operator in the resolvent becomes independent of n x . 
Therefore, the necessary condition for the Laplacian to be universal with Re T is that k z 
dominates k x even for the smallest Reynolds number Re TjYn i n that is considered. This 
poses the above-mentioned aspect ratio constraint on the universal waves. The magnitude 
of the corresponding singular functions increases with Re x J A because of the unit energy 
constraint. The number of waves per unit area and time increases with Re J . In addition, 
the principal singular values of the waves in S ox increase with Re}/ 2 . The above scales 
result in growth of E uu (y; k x , k Z: c) with Re\\ cf. table H 



3.2.3. The outer-in-y class S oy 

For wave speeds close to the centerline, universality of H requires constant \ x /Re T , 
y, \ z , and U c i — c, such that the aspect ratio constraint X x /X z > \^10Re T / Re T ^ 



is 



satisfied; cf. tableland figure [6(g)] As a result, the time T c \ = \ x /{U c i — c) over which 
the wave convects upstream for one wavelength relative to an observer with speed U c \ 
increases with Re T and the convective frequency uj c \ = 2tt/T c i decreases with Re T . The 
aspect ratio constraint follows from similar arguments to those discussed for the outer- 
in-x class. The magnitude of singular functions is independent of Re T since the singular 
functions scale with outer units in the wall- normal direction. The number of waves per 
unit area and time decreases with Re 2 , since their length and convective time increase 
with Re T . The singular values increase with Re 2 and the overall result is that the energy 
density E uu (y; c) increases with Re 2 \ cf. table [Ij 

The waves in the outer-in-?/ class asymptotically approach the streamwise constant 
fluctuations, i.e. n x = 0, as Re T increases. These infinitely long fluctuations exhibit the 



largest linear transient growth in response to initial perturbations in laminar ( Gustavsson 
Butler fc Farrell|1992UReddy fc Henningson|1993l ) and turbulent ( [Butler fc Farrell 



1991 



1993 



.butler & rarrell lyyz; KecLcLy & ttenningso nll^cip ana turbulent ( |£>utler & Farrell 
del Alamo fc Jimenez|2006 Pujals et a/.|2009 ) flows. In addition, they are the most 



highly amplified by the linear dynamics in laminar (|Farrell fc Ioannou 19936 Bamieh 



fc Dahleh|[2QQT| |Jovanovic fc Bamieh|[2QQ5l ) and turbulent ( |Hwang fc Cossu||2010D flows 



subject to stochastic disturbances. 

The effect of Reynolds number on the streamwise constant fluctuations has been stud- 
ied in laminar flows. For example, Gustavsson ( |1991| ) showed that the peak of linear 
transient growth scales with the square of centerline Reynolds number R c i = U c ih/v. For 
the flow subject to harmonic disturbances, Jovanovic & Bamieh ( 2QQ5| ) showed that the 
singular values of H increase as R 2 cl when the temporal frequency uj linearly decreases 
with R c \. No other scales for the singular values were found since the laminar mean ve- 
locity U/U c i = 2y — y 2 is universal with Reynolds number throughout the channel. Our 
study shows that the singular values in the turbulent flow increase quadratically with 
Re T for the waves with defect speeds, U c i — c < 6, and streamwise wavelengths that 
linearly increase with Re T , i.e. > \^10 X z Re T / Re TjYn [ n ; cf. table |ij 



3.3. Universality of the streamwise energy density 

We compute the streamwise energy density of the rank-1 model with broadband forcing 
and illustrate its universal behavior with Re T . These computations build the basis for 
prediction of the streamwise energy intensity at the technologically relevant values of 
Re T in §[4j Because of the unique scales in the inner, outer-in-x, and outer-in-?/ classes of 
wave parameters, we distinguish the corresponding intervals of wave speeds by expanding 
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the premultiplied energy density into the following three integrals: 

r 16 



/-ID 

EuuilJi K z ) = / E uu (y,ti x ,Hi z ,c) 

J5 



dc- 



U cl -6 



16 



u cl (3.4) 

U c i-6 



Similar expansions can be written for E uu (y,K x ) : E uu (y,K z ) : and E uu (y). Notice that 
in spite of the different behavior of singular values and singular functions, the energy 
density increases with Re\ in all three classes of wave parameters. Figure M shows the 
premultiplied one-dimensional energy densities and the energy intensity confined to each 
class of wave parameters and normalized by Re\ . The same contour levels are used for all 
Reynolds numbers, Re T = 3333 (blue), 10000 (red), and 30000 (black), in figures 7(a) 



7(f) Notice that confining the wavenumbers to <S^, S ox , and S oy yields a universal energy 
density as summarized in tables [3] and [4j 

The inner peak of streamwise energy density in the rank-1 model with broadband 
forcing occurs at y + « 11, A+ « 184, and A+ « 44; see figures 7(a) and 7(d) The location 
of the above wall-normal peak represents an integral effect over all wave parameters in 
Si and corresponds to the critical speed c ~ 8.5. The inner peak is comparable with 
the location, length, and spacing of the most energetic structures associated with the 
near- wall cycle , i.e. y + « 15, A+ « 700 — 1000, and A+ « 100 (see, for example, Hoyas & 
Jimenez||2006 ). The outer-in-x peak takes place at \[yy+ ~ 5.2 — 5.8 (c orre spond ing to 

critical speed c ~ U m ), X x ~ 12 — 16, and \f\ z A+ « 20; see figures 7(b) and 7(e) It has 
the same scaling as the VLSMs and its location is comparable with the most energetic 
VLSMs locat ed at \Jyy* 
and at A T 



et al. 2010a 



3.9 and A^ « 6 in boundary layers (see, for example, Marusic 
12 — 20 in pipes and channels (see, for example, Monty et al 



2009). These agreements are striking considering the simplicity of the rank-1 model 



subject to broadband forcing. Figures 7(c) and 7(f) show that the outer-in-y peak takes 
place at y ~ 0.45 (corresponding to critical defect speed U c \ — c w 2) for X x /Re T ~ 0.1, 
and X z « 2. This peak points to much longer structures relative to the LSM structures 
observed in experiments, i.e. \ x ~ 3 (see, for example in boundary layers, Adrian et al 



2000). The same length scales are reported for channels and pipes (see, for example, 
Guala et al.|2006| |Monty et a£|2009| ). 

The one-dimensional energy densities can be integrated in the remaining wall-parallel 
direction to obtain the streamwise energy intensity E uu (y) for the rank-1 model with 
broadband forcing. The dashed curves in figures |7(g)|7(i) are the energy intensities nor- 
malized by Re\ obtained by confining the wave parameters to S ox , and S oy , respec- 
tively. As expected, the energy intensities are independent of Reynolds number when 
confined to the universal classes of wave parameters. The solid curves are obtained by 
integrating the energy density over all wavenumbers and wave speeds 5 < c < U c i, i.e. 
by confining the wave parameters to S e . These figures highlight the selection of two lo- 
cal peaks by the linear amplification mechanism where the inner and outer-in-?/ peaks 
dominate the outer-in-x peak. The inner peak takes place close to the inner peak of the 
streamwise intensity in real turbulent flows. While the energy intensity of real flows ex- 
hibits outer scales near the center of the channel, there is no strong evidence for presence 
of an outer-in-?/ peak even for high Re T . 

As evident from figure 7(g)| the universal inner waves contribute to more than 96% 
of the total energy intensity for y + < 20 for all Reynolds numbers. On the other hand, 
figure 7(i) shows that the universal outer-in-?/ waves capture a smaller amount of the total 
intensity for y = 0.45 as Re T increases; 95% for Re T = 3333 vs. 86% for Re T = 30000. 
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Figure 7. The premultiplied one-dimensional streamwise energy density Re~ 2 E uu (y, X x ) in 
(a)-(c); and Re~ 2 E uu {y, X z ) in (d)-(e); and the streamwise energy intensity Re~ 2 E uu (y), dashed 
curves in (g)-(i) for the rank-1 model with broadband forcing. The wave parameters are confined 
to S = Si in (a), (d), (g); S = S ox in (b), (e), (h); and S = S oy in (c), (f), (i). The Reynolds 
numbers are Re T = 3333 (blue), Re T = 10000 (red), and Re T = 30000 (black). The contour 
levels decrease by 0.05 from their maximum value of 0.25 (a); 0.15 (b); and 0.25 (c), and by 0.1 
from their maximum value of 0.5 (d); 0.4 (e); and 0.8 (f). The solid curves in (g)-(i) are obtained 
by confining the wave parameters to S = S e . 



This is because the aspect ratio constraint in S oy excludes more wavenumbers from S e 
as Re T increases. The excluded waves are not universal with Re T and their contribution 
to the energy intensity is not completely negligible. A similar reasoning explains why the 
universal outer-in-x scale captures 82% of the total energy intensity at \/yy + = VT0 for 
Re T = 3333 vs. 72% for Re T = 30000; cf. figure |7(h)l 

At the end of this section, we recall that the streamwise energy densities and intensities 
thus far were obtained for the model with broadband forcing in A^, \ Z1 and c. In §[4j we 
consider a non-broadband forcing by introducing an optimally shaped energy density. 
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4. Predicting the streamwise energy intensity 

In this section, we introduce a model for predicting the energy intensity of real tur- 
bulent flows by considering a non-broadband forcing in wave speed. This is done by 
incorporating a positive weight function W(c) that amplifies or attenuates the energy 
density E uu (y,c) of the rank-1 model with broadband forcing. Even though W(c) differs 
from a true forcing spectrum (that also depends on the wall-parallel wavelengths), it 
provides the model with sufficient degrees of freedom for predicting the energy intensity. 
In addition, since each wave speed is associated with a certain class of wavelengths, W(c) 
affects different classes of wavelengths as the wave speed changes. 

First, we show that W(c) can be optimally shaped such that the model-based stream- 
wise energy intensity, 

rUai 

E uu ,w(y) = J W(c)E uu (y,c)dc, (4.1) 

matches the intensity of real flows at low Reynolds numbers. Then, we derive similarity 
laws to approximate the optimal weight functions at high values of Re T . These weight 
functions in conjunction with the energy density of the rank-1 model with broadband 
forcing enable prediction of the streamwise energy intensity at technologically relevant 
Reynolds numbers. 

4.1. Optimal weights for small Reynolds numbers 
The weight function W(c) is determined by minimizing the deviation between E UUj w(y) 



in (4.1) and the streamwise energy intensity obtained from DNS, ^w,DNs(y), in the in- 
terval y + > 1 and y < 0.8. We do not enforce matching for y > 0.8 since it requires 
significantly large values of W(c) for wave speeds close to XJ C \. This is because E uu (y,c) 
is considerably smaller and more localized near the centerline compared to other loca- 



tions and results in sensitivity of W(c); see, for example, figure 5(a) for y + > 24000 
corresponding to y > 0.8 for Re T = 30000. Note that the main amplification mechanisms 
for waves with speeds close to U c i is criticality of the singular functions since the non- 
normality effect is small as the mean shear approaches zero. This results in small gains 
and singular functions that are localized in the wall-normal direction. 

Since E uu (y,c) scales with Re 2 T (cf. table H while E uu ^s(y) does not, we find the 
normalized weight function W(c) = Re\ W(c) that solves the following optimization 
problem 

minimize: * nr. 7 s\k + ^ W c ) ™' 

\\Euu,DNs{y)\\i (4.2) 

subject to: W(c) > 0, 2<c<U d . 
Here, ||#(2/)|| e is defined as (note integration in log?/ + ) 

Aog(0.8Re T ) 

\\g(y)\\l = / ^(log^dlog^, 

JO 

to equally penalize the deviation of energy intensities near the inner peak as well as in 
the channel core. The second term in the objective function, 

provides the weight function with smoothness by penalizing the magnitude of W, and 
7w > controls the importance of smoothness relative to matching the model-based and 
DNS-based energy intensities. 
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(d) 



Figure 8. (a)-(b): The optimal weight functions W as a function of c in (a) and U c i — c in (b) 
for Re T = 934 (orange) and 2003 (green). The tangent lines, L, to W(U c i — c) at U c i — c= 3.9 
for each Re T are shown in black, (c) The model-based streamwise energy intensity E uu ,w{y + ) 
for Re T — 934 (orange) and 2003 (gre en) and the DNS-based intensity E uu ^s(y + ) (black) 
are optimally matched by solving (4.2) for each Re T . The respective curves lie on the top of 
each other, (d) The optimal weight functions W for Re T = 934 (oran ge) a nd 2003 (green) are 
compared with the weight functions obtained using the similarity law (4.4) (black dots). 



The optimization problem (4.2) is solved using CVX, a package for specifying and 
solving convex programs in Matlab flCVX Research, Inc.|2012} |Grant fc Bo yd 2008). We 
find the optimal weights fo r the largest Reynolds numbers that are simulated to date 
usin g DNS, i.e. Re T = 934 ( |del Alamo et a/J|2004[ ) and Re T = 2003 ( |Hoyas fc Jimenez 



2006). Even though these are orders of magnitude smaller than the Reynolds numbers for 
which we predict the energy intensity in § |4.2[ they are free of measurement errors and 
useful for finding the optimal weights. We choose j w = 0.2 to strike a balance between 
matching error \\E uu ^s(y) — Euu,w(y)\\e an d smoothness of W(c). The optimization 
problem is robust with respect to the choice of j w . For example, changing j w by a 
factor of two has negligible effect on the matching error and the optimal weights for 
7 < c < U c i — 2 while slightly modifying W(c) elsewhere. 

Figures 8(a) and 8(b) show the optimal weights as a function of c and U c \ — c for 
Re T — 934 and 2003. These weight functions match ^ nn ,DNS an d E UUj w with a relative 
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error of approximately 0.2%; see figure 8(c) As expected, W is qualitatively similar for 
c < 16 and U c i — c < 6 since both the model-based and DNS-based intensities exhibit 
inner and outer scaling in the respective regions. Figure 8(b) shows that W(U c i — c) 
approximately coincides for Re T = 934 and 2003 for U c i — c < 3.9. We denote this 
universal function by Wi(U c i — c). For simplicity, the weights are approximated by linear 
functions in the logarithmic region for 16 < c < U c i — 3.9. These lines are denoted by L 
(black) and intersect for L(c = —2) = 19.88 and L(U c i — c = 3.9) = 2.54. This gives an 
analytical expression for L as a function of wave speed and Reynolds number since U c i 
varies with Re T 

'U d -c-3.r 



L(c; Re T ) = 2.54 + 17.34 



Ud - 1-89 



(4.3) 



As the Reynolds number increases, W(c) is shifted upward for c < 16 by the kick that 
it receives from the logarithmic region. This is expected since the DNS-based energy 
intensity increases with Re r in the inner region while E uu (y + , c) remains constant. More 
discussion about the relationship between the weights in the logarithmic and inner regions 
are provided in § |4.2| Motivated by these observations, we formulate a similarity law for 
the weight function: 



W 2 (c) + (L(c;Re T )-L (c)), 
W(c;Re T ) = I L(c;i?e r ), 



U c i 



2 < c < 16, 

16 < c < U c i ~ 3.9, (4.4) 
3.9 <c< U ch 



that consists of three segments: A universal outer segment represented by W\(U c i — c) 
for U c i — 3.9 < c < Ud] a Reynolds number dependent linear segment L(c; Re T ) that is 
analytically determined by ( |4.3| ); and an inner segment composed of a universal function 
W2(c) for 2 < c < 16 superposed by a linear function L(c; Re T ) — Lq(c) where L (c) = 
L(c;Re T = 934). Figure [8(d)] shows that the optimal weights computed by solving ( |4.2[ ) 
are well-captured by the weights formulated using the similarity law (4.4). 



4.2. Predictions at high Reynolds numbers 



The similarity law in (4.3)-(4.4) is used to predict the weight functions, and consequently, 
the streamwise energy intensity at high Re r using (4.1). Figure [9] shows the predicted 
weights and energy intensities for Re T = 934, 2003, 3333, 10000, and 30000. As evident 
from figures [9(c) | and 9(d) the model-based predictions are consistent with the exper- 
iments of channel flows (Monty 2005) up to Re T = 4000 and boundary layers (Kunkel 
|fc Marus ic 2006) up to Re r = 23013. A logarithmic dependence of the energy intensity 
on the distance from the wall is predicted in the logarithmic region of the mean velocity 
which is consistent with recent experiments ( |Hultmark et a/.||2012 Marusic et a/.||2013 ) 
and predictions of the attached eddy hypothesis (see, for example, Perry fc Chong|1982 ). 
The predicted intensity for Re T = 3333 lies between the experimental data for channel 
flows at Re r = 3165 and 4000 (Monty 2005| ). The experimental measurements are not 
accurate near the wall as they suffer from spatial resolution issues (see, for example, 
Hutchins et al 2009). Figure 9(d) shows that the predicted values of the inner peak are 
consistent with the boundary layer measurements that are corrected (Kunkel & Marusic 



2006) based on the attached eddy hypothesis. On the other hand, our predictions of the 
energy intensity in channel flow are larger than the data in boundary layers close to 
the outer peak in the logarithmic region. This is expected since recent experiments have 
shown that the large structures are more energetic in internal flows such as channels and 
pipes compared to boundary layers (see, for example, Monty et a/.|2009 ). This difference 
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Figure 9. (a)-(b): The optimal weight functions W obtained using the similarity law (4.4) as a 
function of c in (a) and U c i — cm (b). (c)-(d): The corresponding model-based streamwise energy 
intensity E uu , w (y+) for Re T = 934 (orange), 2003 (green), 3333 (blu e), 10000 (red ), and 30000 
(black). The symbols in (c) are experiment al data from channel flows ( |Monty|2005 ) Re T = 3165; 
V and 4000; A, and from boundary layers ( jKunkel Marusic||2006|) lor Re T = 5813; x, 13490; 
+, 23013; The symbols in (d) are the corrected (Kunkel & Marusic 2006) data in the inner 
region using attached eddy hypothesis. The lines show logarithmic scaling. 



was attributed to the observations suggesting that the VLSMs are longer in internal flows 
than in boundary layers. 

Obtaining the results of figure [9] requires computation of the streamwise energy density 
of the rank-1 model with broadband forcing at the respective values of Re r . Alternatively, 
the universal behavior of E uu (y; k x: k z ,c) can be used to avoid these computations. In the 
present study, we employ the universality (invariance with Re T ) of E uu (y\ k x , k z , c) for Si 
to predict the inner peak of the streamwise intensity at arbitrary high Re T . Expanding 
the weighted energy density according to the wave speed and substituting for the weight 
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Figure 10. (a) The red curve is the model-based energy intensit y E uu ,w(y) for Re T = 10 . 
The blue curve is the contribution from the universal function in (4.5), and the black curve is 
the contribution from the inner class of wave parameters Si . (b) The contribution to the energy 
intensity from the inner class of wave parameters for Re T = 10 4 , 10 5 , and 10 6 . 



function using the similarity law (4.4) yields 



E 



uu,W 



(y) 



rl6 
J2 



W 2 {c)(Re- 2 E uu (y,c)) 1 



L(c; Re T ) - L (c) 
W 2 (c) 



dc- 



universal 



rU cl -3.9 rU cl 

/ L(c,Re T ) (Re- 2 E uu (y,c)) tic + / W^c) (Re~ 2 E, 

J 16 JU c i-3.9 



L (y,c)) dc. 
(4.5) 

The first integral, corresponding to the inner class of wave parameters contains a 
universal function multiplied by a coefficient L(c;Re T ) that also appears in the second 
integral for the faster and larger waves in the logarithmic region. It represents the contri- 
bution from the inner class of wave parameters that are coupled with and amplified by 
the large scales in the logarithmic region. This is similar to the model that Marusic et al. 



(20106) proposed to capture the influence of the large scales ul (close to the geometric 
mean of the logarithmic region) on the small scales us close to the inner peak of the 
energy intensity 

us = u*(l J r j3u L ) + au L . (4.6) 

For the purpose of the present study, ( |4.6| ) implies that the small structures are deter- 
mined by a universal inner-scaled function u* multiplied by a coefficient 1 + (3ul that 



increases with the energy of the large structures. Physically, the first term in (4.6) de- 
scribes the amplitude modulation of small scales by the large scales and the second term 
represents the direct superimposition of the large scales on the inner-scaled near-wall 
peak flMarusic etaJ\2QlQb ). 



The blue curve in figure 10 (a) | shows the contribution of the universal function in (4.5) 
to the energy intensity. This is equal to the contribution of the inner class of wave 
parameters to the energy intensity for Re r = 934, i.e. for L = Lq. In other words, the 
inner class of wave parameters is not influenced by the large scales in the logarithmic 
region for Re T = 934. This is expected since Re T « 1000 is the smallest Reynolds number 
where the purely inner and outer scales are separated by a logarithmic region. The black 
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Figure 11. Variation of the inner (y + — 15) and outer (y + — 3.9Re ^ 2 ) peaks of the stream wise 
energy intensity with Reynolds number. The figure is adapted from |Marusic et aTj ( |2010a] ). The 
black symbols are experimental and simulation data from channels and boundary layers. The red 
symbols are the predicted inner peak intensities obtained from the present model for turbulent 
channels. 



curve shows the contribution of the first integral in (4.5) to the streamwise intensity for 
Re r = 10 4 . Notice that the large scales from the logarithmic region increase the inner 
peak by amplifying the universal function through the coefficient 1+(L— Lo)/W2. The red 
curve is the total intensity obtained by integrating the contribution of all wave parameters 
S e . The inner peak is captured by the first integral and the direct superimposition of 
the large scales on the inner peak is negligible. Therefore, the first integral readily yields 
the behavior of the streamwise intensity near the inner peak. For example, figure [lQ(b)| 
illustrates how the more energetic large scales at Re T = 10 5 and 10 6 further increase the 
inner peak relative to Re T — 10 4 by amplifying the universal function 

The red circles in figure [TT] show the predicted value of the inner peak of the streamwise 
intensity as a function of Re T up to 10 10 . These predictions are made at no additional 
cost using the universal energy density for the inner class of wave parameters and the 
similarity law for the weight functions. This figure is adapted from figure 8 in Marusic 



et al. (2010a) where the DNS and experimental data from channels and boundary layers 
are summarized (open and filled black symbols). The black filled squares and circles, 
respectively, show the magnitude of the inner (y + = 15) and outer (y + = 3.9Re^ 2 ) 
peaks in recent boundary layer experiments (Ma rusic et a/.|2Q10a ). Using these data, the 
authors proposed two possibilities for the behavior of the inner peak at high Reynolds 
numbers. The first possibility is to extrapolate following the trend suggested by the 
filled black squares (line 1). The second possibility, motivated by the fact that the large 
scales increase the energy of the small scales, is to extrapolate following line 3 which is 
parallel to line 2 that captures the variation of the outer peak with Re T . Notice that 
the available data is not sufficient to identify the correct behavior of the inner peak 
as Re T increases. On the other hand, the red circles represent a viable model-based 
prediction of the inner peak at high Reynolds numbers. These results are obtained for 
channels and are potentially different than boundary layers. Our ongoing research, to be 
reported elsewhere, shows that in addition to the inner peak, one can also predict the 
streamwise intensity in the logarithmic region at arbitrary high Reynolds numbers using 
the universality of E uu (y; k x ,k z ,c). 
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5. Concluding remarks 

Starting from the NSE, we highlighted the low-rank nature of the resolvent, formulated 
for individual wall-parallel wavenumbers and frequencies, and illustrated its power by 
showing that the most energetic motions of real turbulent flows correspond to wavenum- 
bers and frequencies whose resolvent is approximately rank-1 (in the wall- normal di- 
rection). Motivated by this observation, we studied the streamwise energy density of 
the rank-1 model subject to forcings in the wall-parallel directions and time that were 
broadband and optimized, or "trained" , with respect to the available DNS data. 

It was shown that the resolvent admits three classes of wavenumbers and wave speeds 
where the corresponding principal singular values and singular functions exhibit universal 
behavior with Reynolds number. These classes are directly related to the universal regions 
of the turbulent mean velocity (which is assumed known a priori) and thus are primarily 
distinguished by the wave speed: (i) a truly inner-scaled class of waves with constant 
speeds in the inner region of the turbulent mean velocity; (ii) a class of waves with outer- 
scaled length and constant defect speeds relative to geometric mean of the logarithmic 
region of the turbulent mean velocity; and (iii) a class of waves with outer-scaled height 
and width and constant defect speeds relative to the centerline. 

The proposed analysis effectively divides the streamwise energy density of the rank- 
1 model with broadband forcing into three universal regions and enables its scaling 
to arbitrary large Reynolds numbers. In addition, it was shown that the most energetic 
wave parameters and the corresponding scales roughly agree with the dominant near- wall 
motions in real turbulent flows. A weighted energy density was proposed to account for 
a non-broadband forcing in the wave speed. We illustrated that the weight functions can 
be optimally shaped such that the weighed energy density of the rank-1 model matches 
the energy intensity of real turbulent flows at low Reynolds numbers. 

As the Reynolds number increases, the optimal weight functions increase for wave 
speeds in the inner region of the mean velocity. Representation of the optimal weights 
using similarity laws revealed that the amount of upward shift is linearly correlated with 
the weight function for wave speeds in the logarithmic region of the mean velocity. In 
other words, the weight function increases with the energy intensity of the large scales 
and amplifies the universal inner-scaled energy density of the rank-1 model. Therefore, 
it implicitly captures the well-known coupling of small scales with the large scales and 
their subsequent amplification in real turbulent flows. One of the main results of this 
study is that the rank-1 approximation, together with the optimal weight functions and 
the (well-known) mean velocity profile, is sufficient for predicting the streamwise energy 
intensity at high Reynolds numbers. 

Our ongoing research, to be reported elsewhere, is focused on the analytical expression 
of the streamwise energy density for wave speeds in the logarithmic region of the mean 
velocity and a priori derivation of the weight functions. Addressing the limitations and 
implications of the low-rank model for predicting the wall-normal and spanwise energy 
spectra as well as the Reynolds stresses is another topic of future research. 
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Appendix A. Derivation of the Reynolds number scalings 

We derive the classes of wave parameters that induce universal behavior on the low- 
rank approximation of H. 

A.l. The intrinsic inner class 

Consider the inner coordinate where all spatial coordinates are scaled with the viscous 
unit vju T 

x + = Re T x, y + — Re T y, z + = Re T z. 
The differential operators in y and the wavenumber symbols in the inner coordinates are 

d/d y + = Re^d/dy, k+ = Re' 1 ^, k+ = Re' 1 ^, A+ = Re~ 2 A. 



Consequently in the inner coordinates, the operators Ra, C, and in (3.1) scale as 





Ra 



C 



Re T X 1 
Re 2 T X 3 Re T X 2 

d Re~ 1 C 2 
C 3 
O4 Re r O5 



Y 3 
Re T Cj2 



Re~ Y2 



ct 



cl 



Re T C, 



t 

5 J 



(Al) 



For given n+ and the operators C\ to C5 and their adjoints are independent of 
Re T . On the other hand, the operators X\ to X3 and Y\ to Y3 contain spatially- varying 
coefficients, U — c and its first two derivatives, that depend on Re T . As discussed in the 
beginning of § [3J U scales with y + and is independent of Re T for y + < 100. Therefore, 
for given and c < U(y + = 100) = 16, the operators X\ to X3 and Y\ to Y3 

are independent of i?e r when acting on functions whose supports are inside the interval 
< 2/ + i$ 100- Since the principal singular functions of are localized around the 
critical layer (i.e. the wall- normal location where the turbulent mean velocity equals c), 
the singular functions of H are negligible outside y + < 100 for c < 16 and all of the 



aforementioned operators are effectively independent of Re T . It follows from (Al) that 

Re T 1 Hi2 Re r 1 His 
CR A C^ = Re~ 1 H 2l Re~ 1 H 22 Re~ 1 H 2 s 

Re; 1 H 32 Re^Hss 



H 



Re~ 1 Hu 

Re~ 1 H 2 i 
Re^H 31 



where the operators are effectively independent of Re T when acting on their principal 
singular functions. Therefore, the princip al sin gular value of H is proportional to Re~ x . In 
addition, the orthonormality constraints (2.9) on ip 1 and (j> 1 require that these functions 
scale as Re}/ 2 . This is because the supports of ^ 1 and cj) 1 are independent of Re T in 
inner units (hence, proportional to Re~ x in outer units). In other words, ^i(y) and 
<j>i(y) become thinner and taller as Re r increases. Finally, the streamwise energy density 
Euu = k 2 k z o~ 2 \ui\ 2 scales with 

(Re r ) 2 (Rer) (Re' 1 ) 2 \Re\' 2 ) 2 = Re\. 

A. 2. The intrinsic outer-in-x class 

Consider the outer-in-x coordinates where x is scaled with the outer unit h and y and z are 
scaled with the geometric mean of the inner and outer units. The wall-normal differential 
operators and the spanwise wavenumber symbol in the outer-in-x coordinates are given 

by 

d/dy = Re~ 1/2 d/d y , k z = Re~ l/2 n z . 
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For given k x and k z , the Laplacian 

A = Re T dyy 



Re T (k z ) 2 , 



approximately scales with Re T if Re T {h z ) 2 dominates k, 2 for all values of Re T . It suf- 
fices that this condition is satisfied for the smallest Reynolds number Re T ^ n that is 
considered. This requires 



X x 



K z 



> 



10 



In the outer-in-x coordinates, the operators Ra, C, and scale as 

-l 

X , SI 

Ra 



(A 2) 



C 



X Y 
Re T X 3 X2 

Re- 1/2 Ci Re~ 1/2 C 2 
C 3 _ 

O4 Re r C5 



Y 1 
Re T % % 



Re- 1/2 C\ C\ C\ 
Re\ /2 C\ C\ 



(A3) 



For given k x and k z that satisfy the constraint (A 2), the operators C\ to C5 and their 
adjoints are approximately independent of Re T . We recall from figure [6(c)] that U (y) — U m , 
with Um denoting the mean velocity at the geometric mean of the logarithmic region, is 
independent of Re T for 50/ \/Re T ^ n < y < 0.2y/Re T ^ m i n . Therefore, for given k x , k z: 
and \U m — c| < (l/2ft) log(i?e Tjm i n (50/0.2)), the operators X\ to X3 and Y\ to I3 are 
approximately independent of Re T when acting on functi ons w hose supports are inside 
the interval 50/ \jRe T ^ n < y < 0.2^i?e Tjm i n . Following (A3) and using the locality of 



the principal singular functions of H around the critical layer, we have 



H = CR A C ] 



Hu 

Re~ 1/2 H 21 

Re; 1/2 H 3 i 



Re x J 2 H 12 Rel /2 H 13 



H22 
H32 



H23 
H33 



where the operators Hij are approximately independent of Re T when acting on their 
principal singular functions. Therefore, the principal singular value of H is pro port ional to 
Re\l 2 In addition, since and <j) 1 scale in £/, the or t honor mality constraints (2.9) require 
that these functions be proportional to Re]J A . Then, the streamwise energy density E uu = 
K 2 K z a 2 \ui\ 2 scales with 



Ret 



A. 3. The intrinsic outer-in-y class 

Consider the outer-in-?/ coordinates where y and z are scaled with h and x is scaled 
hRe T . The streamwise wavenumber symbol in the outer-in-?/ coordinates is given by 
k~ = Re T K x . The Laplacian 

A = dyy — Re T (k x ) 2 — k z , 

is independent of Re T if k 2 dominates Re~ 2 (n~) 2 for all values of Re T . For fixed n~ and 
k z , it suffices that 



tv x Re T 



10 



(A4) 
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C 








x 3 


Re- 1 ^ _ 


Re~ 1 C l 


c 2 


c 3 





c 4 


Re~ l C z _ 



i?e T Yi 

Re-'Gl Ct 
C\ 



c2 

ReZ l Cl 



(A 5) 



For given and n z that satisfy the constraint (A 4), the operators C\ to C5 and their 
adjoints are approximately independent of Re T . In addition, the defect velocity U c i — U(y) 
is independent of Re T for y > 0.1. Therefore, for given ac~, and U c \ — c < C/ c / — U(y = 
0.1) = 6, the operators Xi to X3 and Yi to I3 are approximately independent of i?e r 



= CR A C ] 



Re 2 T H 12 
Re T H 22 
Re T Hs2 



Re 2 T Hi3 
Re T H23 
Re T Hs3 



(A 6) 



when acting on functions whose supports are inside the interval 0.1 < y < 1. From (A 5), 
we have 

^31 

Owing to the locality of the principal singular functions around the critical layer, the 
operators are approximately independent of Re T when acting on their principal 
singular functions. Therefore, the principal singular value of H is proportional to Rej . 
Since and <j) 1 scale in the outer length scale, the orthonormality constraints (2.9) 
require that these functions be independent of Re T . Finally, the streamwise energy density 
Euu = KxKz°~i\ u i\ 2 scales with 



{Re-y(l) (RelfdY = Re 
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